On the Kirzhnits gradient expansion in two dimensions 
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We derive the semiclassical Kirzhnits expansion of the D-dimensional one-particle density matrix 
up to the second order in h. We focus on the two-dimensional (2D) case and show that all the gradient 
corrections both to the 2D one-particle density and to the kinetic energy density vanish. However, 
the 2D Kirzhnits expansion satisfies the consistency criterion of Gross and Proetto [J. Chem. Theory 
Comput. 5, 844 (2009)] for the functional derivatives of the density and the noninteracting kinetic 
energy with respect to the Kohn-Sham potential. Finally we show that the gradient correction to 
the exchange energy diverges in agreement with the previous linear-response study. 
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I. INTRODUCTION 

Gradient expansions provide a natural path to cor- 
rect the local-density approximation (LDA) for slowly- 
varying densities as already suggested by Hohenberg and 
Kohn in their seminal paper on density-functional the- 
oryi (DFT). The second-order gradient expansion for the 
kinetic energy was shown to be very useful, but on the 
other hand, systematic gradient expansions for the ex- 
change and particularly for the correlation energies faced 
problems that were later corrected - at least from the 
practical viewpoint - by generalized-gradient approxima- 
tions-^ (GGAs). Nevertheless, gradient expansions pose 
still open questions, especially in reduced dimensions 
such as the two-dimensional electron gaa^ (2DEG). The 
interest in the 2DEG arises from a multitude of applica- 
tions in, e.g., quantum Hall and semiconductor physics. 

Semiclassical gradient expansions can be regarded as 
alternatives to the standard approaches based on Tay- 
lor expansions and linear-response formalism. Although 
semiclassical methods do not give access to the correla- 
tion energy, they can be used to derive simple density 
functionals for the Kohn-Sham (KS) kinetic energy Tg 
and the exchange energy density (Ref. 0). Here we 
focus on the semiclassical Kirzhnits commutator formal- 
ism;^ It has been previously used to derive the lowest- 
order (second order in h) gradient correction terms to 
the one-particle density matrix 7(r, r') and to ex in three 
dimensions (3D))^ as well as for Tg in D dimensionsi^ 
Higher-order corrections in 3D have been considered in 
Ref. 0. Salasnicbi found that in 2D the gradient correc- 
tions to the kinetic energy vanish, which was in agree- 
ment with earlier results based on the response-function 
approach/^ This 2D gradient correction was also recently 
studied for systems at finite temperature. 

In this paper we use the Kirzhnits method to derive 
the lowest-order gradient corrections to the one-particle 
density matrix in D dimensions. Then we focus on the 2D 
case and show that all the corrections to the one-particle 
density n(r) vanish, and, in agreement with Ref. they 



vanish also for Tg . Due to the resulting simple expressions 
for n(r) and Tg the consistency criterion of Gross and 
Proettcii that couples the functional derivatives of Tg 
and n{r) is trivially satisfied. Finally, we show that the 
gradient corrections to ex diverge in the 2D Kirtzhnits 
expansion, which is in agreement with the linear-response 
results of Gumbs and Geldart 



II. KIRZHNITS EXPANSION IN D 
DIMENSIONS 

The exchange energy Ex and the KS kinetic energy Tg 
can be expressed 




where the one-particle density matrix can be written in 
terms of the Fermi energy e^:^ as 

7(r,r')= 

= ^e(eF - ej)ip*{r)ipj{r') 
j 

^{r\e{e^-i^vs)W). (3) 

Here ipi are the solutions of the single-particle KS equa- 
tion, t is the kinetic energy operator, Vg is the KS po- 
tential, and is the Heaviside step function. Now we 
define the local Fermi energy Ep = ep ~ Vg{r) and use 
the plane- wave decomposition as 

7(r,r') = E / d''k{r\Q{EF-i)\ka){ka\r'), (4) 

where |ka) (with a as the spin index) are eigenfunctions 
of the momentum operator p. 



We introduce the abbreviated notations: <d{EF "i) = 
f(a + b), f = 0, a = -i = -f/2, and b = Ep = 
kp/2. Now we can use the inverse Laplace transform, 
the Fourier-MeUin integral, to show that the operator 
Q{Ef — t) acts on eigenfunctions |k) as 



expression of the density matrix. 
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where c = Re(/3) > is arbitrary, but chosen such that 
the contour path of the integration is in the region of 
convergence of F{(3). The commutation problem of op- 
erators a and b can be avoided by introducing a new 
operator K{I3) (Refs. and[^ such that 



t3a 



= eP'K{f3)e 
Thus, we obtain for Eq. ([5]) an expression 
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where the operator a has now been replaced with eigen- 
value a in the exponential function, so that it commutes 
with the operator b. 

The expression for the operator K{(i) is obtained by 
expanding it in a power series with respect to /3, 



(8) 



n=0 



Differentiating both sides of Eq. ^ with respect to /3, 
and expanding all exponential functions in a Taylor se- 
ries, leads to the recurrence relatior^ 
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Inserting Eq. ^ in Eq. ^ yields 



f{a + b)\a)^Y. 
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where /^"^ is the n:th derivative of the function /. The 
Kirzhnits expansion in Eq. leads to the following 



oo ^ 

^ — n ^ \- ^ 



n—2 Q— ± 



Ef 



(r|d„|ka)(k«|r'). 



(13) 



where 5^^^ is the n:th derivative of the delta function. 
The first-order term of Eq. ([T^ 7^°^ corresponds to the 
zeroth order solution of the one-particle density matrix 
and it can be written as 
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where dfl is the {D — l)-dimensional angular volume el- 
ement and is the angle between vectors y and k. This 
term generates the exact exchange energy for the ho- 
mogeneous electron gas which can be used as the LDA 
in an inhomogeneous system. Note that in Eq. ([T4| we 
have defined the relative and center-of-mass coordinates 
as y = r — r' and R = (r -|- r')/2, respectively. 

Higher-order terms of the Kirzhnits expansion 7^"^ can 
be determined by calculating higher derivatives of the 
delta function and multiple commutators of Ep and t 
that lead to multiple derivatives oi kp. The second-order 
(V^) inhomogeneity correction consists of three terms, 
ry(^)[r, r') = B2 + B3 + B4, where 
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Here we use a definition z = z(R, ?/) = fci?(R)|y| and 
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expressions 



III. TWO-DIMENSIONAL CASE 





= J d^kS' 










},d.-4 


2)I{z) + zI'{z)]; 


(18) 




^ J d^kS" 


Ef - 

2 


giky 





L.d-6 



[d^ - 6d + 8] I{z) + [2d - b)zl\z) 



+ z^I"{z) 
h{z) = j d^kS'" 



(19) 



Ef - -k' 



o«ky 



kp 



16 



[d^ - 12<i2 + Ud - 48] /(z) 
+ 3{d^ -7d+ n)zl\z) + 3{d - 3)z^I"{z) 
+ z3/(3)(z)l (20) 



Combining our results leads to the semiclassical expan- 
sion of the density matrix of the form 



7(r,r') =7^°nr,r')+7^'Hr,r') 



- S„,„, i A + BiVRk'p) • ^ + C V^fc^ 
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where A, B, C, D, E, and F are given by 
2 
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Using these equations it is straightforward to proceed 
with the calculation of the one-particle density matrix 
and the exchange energy density in 2D. 



A. One-particle density matrix and the kinetic 
energy 

In the two-dimensional case we obtain the following 
expression for the remaining integral in the one-particle 
density matrix in Eq. (1141) . 



I{k) = / dl]e*'=^'=°'^^ = 27rJo(z). 
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This leads to 
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where Jn{z) is the Bessel function of the first kind in 
order n. This expression is the central result of this paper 
and in the following it is used for further analysis. 

We first notice that the one-particle density has a sim- 
ple form 



Ztt 



(24) 



i.e., all the gradient corrections vanish. This may seem 
as an unexpected result in view of the known gradient 
expression in 3D.^ Secondly, we calculate the KS kinetic 
energy by inserting Eqs. ((23|) and (p4|) into Eq. ([2]) and 
find 

T,^ J drt,(r) 

= -|;/<'r{vh.(r,OU, 

/ <at I fiva + Vy) T(R.!,)l , (25) 

I J y=0 



2m 



where is the noninteracting kinetic energy density. Af- 
ter lengthy but straightforward calculations we find 



?52 

ts{r) = Tj-T^n (r), 
2m 



(26) 



which is equal to the Thomas-Fermi expression. Again, 
the gradient correction (von Weizsacker term) is zero. 
This 2D result is in agreement with previous Kirzhnits 
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expansion for the _D-dimensional kinetic energj«i as well 
as with results obtained using alternative methods i^ii^ii^ 
In their recent work Proetto and Gross^i have derived 
a rigorous condition to test the consistency of approxi- 
mations made for the density and the KS kinetic energy. 
The condition is given by 



dr'v^ir')'4^, (27) 



Svs{r) J Svs{r) 

where Vg is the KS potential. We note that the condi- 
tion follows from the Euler equation minimizing the KS 
energy, i.e., 

+ (28) 

where n is the chemical potential. Multiplying both sides 
with 5n{r')/Svs{r) and integrating over r' directly yields 
Eq. ([271). The condition means that STs[n\/Sn{r') = ep- 
Vg must be also valid. Using Eqs. ((M)) and ([^S)) . and 

kp = \J1m[ep — Ws)/fi2 we find 

2m 



— — Trrt — - — h^p = ep ~ Vg 



(29) 



6n(r') 

Thus, Eq. dill) is fulfilled for the 2D (and also 3D) resuhs 
of the semiclassical Kirzhnits expansion. 



B. Exchange energy 



In other words, the exchange energy density is clearly di- 
vergent in the 2D Kirzhnits expansion. Our result agrees 
with the finding of Gumbs and Geldartj^^ who used per- 
turbation theory and linear-response formalism to derive 
the second-order gradient terms for both the kinetic and 
exchange energies in D dimensions. They arrived at the 
same result also by using the Wigner-Kirkwood expan- 
sion. Hence, as confirmed in this work from the semi- 
classical point of view, the divergence of the systematic 
gradient expansion for the exchange energy seems to be 
an inevitable mathematical fact. However, to the best of 
our knowledge, the underlying physical reason that makes 
the 2D situation specially divergent in contrast with the 
ID and 3D cases remains unknown. We hope that the 
present analysis encourages further examinations from 
that viewpoint. 

The divergence of the exchange energy in 2D can 
be considered unfortunate in view of functional devel- 
opments in 2D, although first GGAs in 2D have al- 
ready been obtained,— and several other 2D function- 
als have been derived, for example, in the framework 
of meta-GGAsJ^ A natural next step, as already dis- 
cussed in Ref. [l^, would be considering expansions in 
quasi-2DEG by introducing a finite width of the system. 
This would resemble also the experimental situation in 
low-dimensional nanostructures such as in semiconduc- 
tor quantum dots. 



Knowledge of the gradient corrections to the one- 
particle density matrix in Eq. (|23p immediately moti- 
vates to search for an expression for the exchange energy 
defined in Eq. ([1]). Using Green's first theorem we ob- 
tain the second-order expansion of the exchange energy 
density in H in terms of the gradients of fci^: 



1927r k% 



dz D{z), 



(30) 



where 

D{z) = z2jo(z)2 + (z2 _ 4) Ji(z)2 - 2zJi{z)J2{z) 

+ 4zJo(z)Ji(z) + z2jo(z)J2(z). (31) 

Expanding Bessel functions in Taylor series and using a 
regularization of divergent Coulomb integrals leads to 
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IV. SUMMARY 

In summary, we have derived the second-order gra- 
dient corrections to the one-particle density matrix in 
the semiclassical Kirzhnits expansion in D dimensions. 
In two dimensions the corrections vanish in the diago- 
nal of the density matrix, i.e., in the one-particle den- 
sity. Similar vanishing occurs in the noninteracting ki- 
netic energy in agreement with Ref. 0, and leads to the 
fulfillment of the consistency criterion of Ref. [ll|. Fi- 
nally, we have shown that the exchange energy of the 
two-dimensional Kirzhnits expansion diverges in agree- 
ment with the linear-response theory. We hope that the 
present work motivates further attempts in the system- 
atic derivation of gradient corrections in the quasi-two- 
dimensional electron gas. 
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